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Abstract 

We will study entangled two-photon states generated from a two-mode supersymmet- 
ric model and quantify degree of entanglement in terms of the entropy of entanglement. 
The influences of the nonlinearity on the degree of entanglement is also examined, and 
it is shown that amount of entanglement increase with increasing the nonlinear coupling 
constant. 
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1 Introduction 

Perhaps, quantum entanglement is the most non-classical features of quantum mechanics 
which has recently attracted much attention although it was discovered many decades ago 
by Einstein, Podolsky, Rosen [1] and SchrOdinger [2]. It plays a central role in quantum 
information theory and provides potential resources for communication and information pro- 
cessing [3, 4, 5]. By definition, a pure quantum state of two or more subsystems is said to 
be entangled if it is not a product of states of each components. A lot of works have been 
devoted to the preparation and measurement of entangled states. Moreover the possibility for 
generation of the entangled states with a fixed photon number has been theoretically studied 
[6, 7, 8, 9]. Duan et al described an entanglement purification protocol which generates maxi- 
mally entangled states with fixed photon number from squeezed vacuum states or from mixed 
Gaussian continuous states by the quantum non-demolition measurement [6, 7]. Quantum 
teleportation using an entangled source of fixed photon number has also been theoretically 
investigated in [8] . Liu et al are used a system of two coupled microcrytallites as a source with 
fixed exciton number and quantified entanglement of the excitonic states [9]. Therefore, the 
generation of a new entangled source with fixed photon number is an interesting task both 
from experimental and theoretical viewpoints. 

In this contribution, it is shown that a two-mode field with a two-photon interaction can 
be used as a good source for generation of entangled states with fixed photon number. We 
will study entangled states generated from two degenerate bosonic systems with fixed photon 
number, and we concern on quadratic nonlinearity between modes to use Higgs algebra as 
the spectrum generating algebra of the corresponding Hamiltonian [10, 11]. We also restrict 
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ourselves to the case that total number of photons is odd. For this case, Debergh in [10] have 
shown that the corresponding Hamiltonian is supersymmetric [12]. 

A number of entanglement measures have been discussed in the literature, such as the 
von Neumann reduced entropy, the relative entropy of entanglement [13] and the so called 
entanglement of formation [5]. In order to discus entanglement of the states, we use von 
Neumann reduced entropy which has widely been accepted as an entanglement measure for 
pure bipartite states. 

The organization of the paper is as follows. In section 2 we introduce a quantum optics 
model for two bosonic system with supersymmetric feature. An analyticl solution of the 
Hamiltonian is also given following the method of Ref. [10]. In section 3 the analytical results 
of section 2 are employed to generate entangled two-photon states with fixed photon number. 
Some examples are also considered in section 3. The paper is concluded in section 4 with a 
brief conclusion. 



2 The two-mode supersymmetric Hamiltonian 

In this section we shall introduce and analyse a model for nonlinear interaction between two- 
mode field. Our method is based on the analysis given in Ref. [10]. Let us consider the 
following family of Karrassiov-Klimov Hamiltonian [14] which describes multi-photon process 
of scattering, i.e. 

H = u)ia\ai -\- 0:202,0.2 + fi'(«i)*«2 + 9*o\{a2)'^ (1) 

where < r < s, 5 is coupling constant and Wi {i = 1, 2) refer to angular frequencies of two- 
mode field characterized by annihilation and creation operators Ui, a\ respectively, satisfying 
[aj,a]] = 5ij. Hamiltonian (1) can be rewritten as 

H = {uji-\- u)2)Ro + (swi - ruj2)Jo + gJ+ +g*J-, (2) 

where [11] 



and 



-Ro = — - — (ra\ai -\- sala2), (3) 
r -\- s 

Jo = (ojai - 4a2), J+ = {a\ya^, J- = af (4)''. (4) 
r -\- s 



It can be easily show that 
and 



[Ro,Jo] = [Ro,J±] = 0, (5) 
[Jo,^±]=±^±> (6) 



for arbitrary values of r and s. 

It is obvious that Ro is a constant of motion, and the total photon number of the two-mode 

system is conserved. Moreover, the infinite dimensional vectors {|ni,n2) = 

n'i,n2 = 0, 1,2,---} are eigenvectors of Rq with corresponding eigenvalues j = ''"^.^^"^ . 
Debergh [10] has shown that in order to have Higgs algebra as the spectrum generating 
algebra of the Hamiltonian (2), we have to add to Eq. (6) the following requirement 

[J+,J_] = 2Jo + 8/3j3, (7) 
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and have shown that [10, 11] this is possible only for r = s = 2, with parameter f3 given by 



4j2 + 4j-2' 
These values of (3 lead to the relations [15] 



i = o,-,i,-. 



m , 



Jslj^'m) = -^\j,m), 



(8) 



(9) 



J±\j,m) = V(j T m){j ± m + T m - ±m + 2)\j,m± 2), (10) 

for TO = — j, — j + 1 ■ ■ ■ . +.7. For a fixed total photon number A'', the vectors \ j, to) are related 
to the two-mode Fock states by 



|j,TO) = |TOi)A|m2)s = 



{a\y+"'{aly-'^ 
V^(j +TO)!(j - to)! 



|0)a|0)b, 



(11) 



where \mi) A®\iTi2) B represent Fock state with toi = j+to photons in mode A and TO2 = j—m 
photons in mode B. 

In order to have more symmetry in Hamiltonian (2), let us suppose that coi = CO2 = 
and concern on the case that g is real. In this case Hamiltonian (1) reduce to 



H = Lj{a[ai + 402) + g ((ai)^ai + af (4)^) 
= 2ujRQ + g{J+-i-J-). 



(12) 



Now, by expanding eigenvectors of (12) as \ipk) = Sm=-j^™ using eigenvalue 

equation H lip k) = Ek\^l>k) , get 



EkC^^ =2jujC'^^ 

+gd^U VU + m)iJ +m- - m + l)(j - m + 2) 
+gd^i2 VU - m)U - m - l)(i + to + + to + 2). 



(13) 



Moreover in order to have supersymmetric Hamiltonian, Debergh concerned on the case that 
j is a half-integer, which leads to twofold degeneracy of all eigenenergies as 



Ek = 2uj + gXk, 



k = l,2,--- ,j + 



1 



(14) 



where Afe is anyone of the j + | different solutions of [10] 



[F{Ak,j,X)r^ 



E 



a2 a2 _ a2 a2 



fe=i 



ife<;, \k-l\Tt2 



Y: AlAfAl-Y^AlA^.A 

yk<l<p, |fe-i|^2, \k-p\^2, \l-p\^2 fe=l 

where Ak are defined by 



Ak = {k{k + l){2j-k){2j-k + l))^, 



k = l,2,--- 



(15) 



(16) 
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Let us denote two eigenvectors of H corresponding to twofold degenerate eigenvalue Ei~ with 
1^^^^) and ivf Now, since Eq. (13) relates coefficient C^^ to C^ja and d^l^, we can, 
without lose of generality, write these two orthonormal eigenvectors belonging to eigensub- 
space Sk as 



l^i'^) = EU Cl%n\jJ - 2n), ^ 



l€^> = EnJ Cj^k-ili-i - 2n - 1), C^%„_, 



(17) 



where 



"j ~ "j-2n— AiA3 - A2„_i ' >J 2' 

J-Cs) _ 1 J>(fe) _ FiA2p,n-^,Xk) , . _ 1 

~ ' "j-2n-l — A-iAi-A-ir, ' -^r-'iJ 2' 



(18) 



where function F has been defined in Eq. (15). In two-mode Fock space representation, Eq. 
(17) can be written as 



14'^) = Erio 2n|2i - 2n)A|2n)s, 

IV-f ) = ES Cfl^n\^3 - 2n - l)^|2n + l)s. 

Finally, evolution operator U (t) takes the following form 



(19) 



uit) = J2 + IV'f l) . (20) 



A;=l 



3 Two-photon entanglement 

In this section we will study entangled states generated by Hamiltonian (12). The entangle- 
ment measure that we are going to use is, the so called von Neumann entropy of reduced 
density matrix which has most widely been accepted as an entanglement measure of pure 
state of a bipartite system. Let lip) be a pure state of a bipartite system with state space 
Ha ® Hb- Entanglement of lip) is defined by 

E{m = -TripAlnpA) = -Tr(pBlnpB) = ^ A^lnA^ , (21) 

n 

where pA is reduced density matrix of subsystem A which is obtained by tracing out subsystem 
B, i.e. Pa = TrsdV') (?A|), Pb is defined similarly, and A„ are square root of nonzero eigenvalues 
of pA and Pb- They are also Schmidt number of state \tl>), i.e. 

IV') ='^K\Un)A\Vn)B, (22) 

n 

where {|w-n)} and {l^n)} are orthonormal states of two subsystems A and B, respectively. The 
definition is based on the fact that although entropy of a pure state is zero, but von Neumann 
entropy of each subsystem is zero only when the state Itjj) is a product state. 
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In this paper we shall consider the case that the total number of photons in the whole 
system is fixed by the initial condition N = 2j, and system is initially in product state 

\m) = \N-L)A\L)B, (23) 

where represents initially N — L photons in mode A and L photons in mode B. By taking 
account of Eqs. (19), (20), (23), we obtain, up to an overall phase factor e*^'^, the final state 
of the system by 

Af-l 

\^(N-L,L)^t)) = J2 an\N -2n- AL)A\2n + Al)b, (24) 



n=0 



where coefficients a„ are defined by 



«"=E^-^'^*^SL^SL-A.' (25) 

fe=i 

and Al is difined such that it is zero (one) when L is an even (odd) integer. Obviously, Eq. 
(24) represents final state of the system in Schmidt form and, accordingly, the von Neumann 
entropy of the reduced density matrix can be obtained easily by 

£;(^-^.^)(t)=-^|a„|2ln|a„|2. (26) 

n=0 

Finally, it should be stress that according to Eq. (26) maximal entangled state a of system 
with the total photon number N is 

= E |iV - n)^|n)B, (27) 



n=0 



where in this case entropy of entanglement is equal to e'^m]^-^ = ln(N On the other hand, 
for state given by Eq. (24), maximum entropy of entanglement is obtained when o„ = y'^^. 



I.e. 



IV'K"''^) = J ^ ^ |iV - 2n - Az.)A|2n + Az.)b, (28) 

where we find ^'^^x '^'' = ln( ''^^-'" ). This means that for a system with fixed photon number 
N, maximum entanglement that can be achieved from Hamiltonian (12) is less than maximum 
entanglement that can be obtained from a system that linear interaction between modes is 
also considered. The difference between these two maximum is, of course, constant and equal 
to ln(2). 

In the rest of this section we will consider some examples in N = 1,3,5,9 and discus 
results. 

1. j = §• In this case the total number of photons of system is 1, and because of the twofold 
degeneracy of eigenvalues, the whole state space of system coincide with eigensubspace 
of the only eigenvalue. Accordingly the final state |V'(i)) differs with initial product 
state only in a total phase factor, therefore, we can not have entanglement. 
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2. j = §• In this case the state space of system decomposes into two eigensubspces, with 
eigenvalues 

Ei=Su)-h VUg E2 = 3u)- Vl2g. (29) 

By starting with two initial states \i>{0)) = \3)a\0)b and IV'(O)) = |2)^|1)b we obtain, 
respectively 



{^^(■''•^^t)) =cos{VT2gt)\3)A\0)B-ism{VT2gt)\l)A\2)B, (30) 
|V'(2'i)(t)> =cos(Vl2gi)|2)A|l)s-i sm{Vl2gt)\0)A\3)B. (31) 

By using Eq. (26), we obtain same value for entanglement of the above two states as 

= - cos^ {Vl2gt) ln(cos2 (VT^gt)) - sin^ iVT2gt) ln(sin2 {Vl2gt)). ^ ' 

Equation (32) shows that entanglement has zero value when ^ = ||- (for A; = 0, 1, • • • ) and 

it takes maximum value ln(2) at times t = ^"^^^^^^ (for A: = 0, 1, • • • ). This, obviously, 
shows that the survival time of maximum entanglement decrease with increasing of the 
nonlinear coupling constant g. 

3. j = §• This case corresponds with a system that has five photons and the state space 
of system decomposes into three eigensubspces, with eigenvalues 

El = £2 = 5a; + As/lg E3 = 5uj - AVlg. (33) 

In this case by considering the initial state as anyone of |V-'(0)) = |5)a|0)_b, |'(/'(0)) = 
|4)a|1)b and |V^(0)) = |3)a|2)b, wc find, respectively, the final state of the system as 

IV'(''°HO) = i^(9 + 5cos(4V75i))|5)A|0)B 

-z^Jsin {AV7gt)\3)A\2)B + (-1 + cos {^Vlgt)) |1)a|4)b, ^^^^ 

\i>^'''\t)) =^(5 + 9cos(4V7,gt))|4)^|l)B 

-z^sin(4V7.gt)|2>A|3)s + 2^(-l + cos(4V7(7i))|0)^|5>B, ^ ' 

|^(3,2)(^)) = -i.f^ sin {AV7gt)\5)A\0)B 



+ cos {4.V7gt) 1 3) A 1 2) B - i ^ sin (AVlgt) 1 1) a 1 4) b . 

Figure (1) demonstrates the evolution of the entropy of entanglement as a function of 
gt for three different initial states with different nonlinear coupling constant g. The 
figure is plotted such that the top horizontal line of each curve corresponds to the 
maximum cntangienicnt ln(3). The maximum ontangiomcnt that can be obtained by 
system is different for different initial state and the system can reach, approximately, to 
maximum entanglement ln(3) only in the case that the initial state is \ip{0)) = |5)a|0) b- 
The Fig. (1) also shows that the siu'vival time of maximum entanglement decrease 
when the difference between photon numbers of two modes A and B of the initial state 
is decreased. As the horizontal axis of the curves is product of coupling constant g 
and time t, it is obvious that by increasing the nonlinear constant g, survival time 
of maximum entanglement decreases. Equations (34), (35) and (36) show that if the 
nonlinear coupling constant g is equal to zero, then \ip{t)) = \tp{0)), that is we can not 
have entangled state. 



Two-Photon Entanglement in a Two-Mode Supersymmetric Model 




10 15 20 

gt 



Figure 1: E^^'^\ E'^^'^^ and i?(3,2) g^j,g plotted as a function of gt in interval [0,2] (curves (a), 
(c) and (e)) and in interval [0,22.015] (curves (6), (d) and (/)). 



4. j = f • Finally we consider as the last example the system with nine photons and ac- 
cordingly the state space of system decomposes into five eigensubspces, with eigenvalues 




v/792 - 24\/56Tg 
v/792 - 24\/56Tg. 



(37) 



The evolution of the entropy of entanglement as a function of gt for four different initial 
states |V(0)) = |9)a|0)s, |V(0)) = |8)a|1)b, |V'(0)) = |6)a|3)b and |V(0)) = |5)a|4)b is 
demonstrated in Fig. (2). The maximum entanglement that can be obtained by system 
is different for different initial states (the top horizontal line of each curve corresponds 
to the maximum entanglement ln(5)). We find that the maximum entanglement ln(5) 
is obtained, approximately, only in the case that there are nine photons initially in one 
of the modes (for example mode A), i.e. [■0(0)) = |9)a|0)b- The survival time of the 
maximum entanglement decrease by increasing the nonlinear constant g and it is also 
decrease by decreasing the difference between photon number of two modes A and B of 
the initial state. 
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(a) (b) 




0.5 1 1.5 2 0.5 1 1.5 2 

gt gt 

Figure 2: E''^'°\ E^^'^\ E'-^'^'> and E'-^''^'> are plotted as a function of gt (curves (a), (6), (c) 
and {d)). 

4 Conclusion 

We studied entangled states generated from two-mode supersymmetric model with fixed pho- 
ton number. We found that only in the case that system has = 3 photons, the maxi- 
mum entanglement can be obtained exactly. For other systems with total photon number 
greater than three, we found that the maximum entanglement is obtained, approximately, 
only in the case that all photons are initially in one of the modes, i.e. |V'(0)) — \N)a\0)b or 
IV'(O)) = \0)a\N)b- The influences of the nonlinearity on the degree of entanglement is also 
examined, and is shown that survival time of maximum entanglement decrease by increasing 
the nonlinear coupling constant g. It is also shown that the survival time of maximum entan- 
glement decreases when the difference between photon number of two modes A and B of the 
initial state is decreased. 
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